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Abstract:  The  present  paper  deals  with  a  theoretical  and  experimental  study  of  a  new  generation  of 
water  desalination  unit  by  solar  energy  using  the  humidification  and  dehumidification  (HD)  principle  is 
constructed  at  the  national  engineering  school  of  Sfax  (34N,  10E),  Tunisia.  The  good  quality  of  distilled 
water  obtained  by  this  new  concept  favours  its  use  for  producing  water  for  drinking  and  irrigation.  A 
mathematical  model  based  on  heat  and  mass  transfers  in  each  component  of  the  unit  is  developed.  The 
resulting  ordinary  differential  systems  of  equations  are  transformed  into  a  system  of  algebraic  equations 
using  the  orthogonal  collocation  method  (OCM)  and  simulated  using  C++  software  in  a  steady  state 
regime.  The  numerical  model  is  used  to  investigate  the  thermal  performance  of  this  kind  of  installation 
exposed  to  a  variation  of  the  control  parameters.  The  thermal  performance  was  evaluated  by  the  gained 
output  ratio  (GOR)  and  the  efficiency  of  the  water  solar  collector.  A  series  of  experiments  was  conducted 
and  compared  with  the  simulation  results  to  validate  the  developed  models.  As  a  result,  the  proposed 
models  can  be  used  for  sizing  and  testing  the  behaviour  of  such  a  type  of  desalination  unit. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

In  many  places  world  wide  drinkable  water  is  already  a  scarce 
good  and  its  lack  will  rise  dramatically  in  the  future.  Missing 
energy  sources  and  no  grid  connections  complicates  the  use  of 
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standard  desalination  techniques  in  these  places.  Today,  sea  and 
brackish  water  desalination  plants  are  well  developed  in  industrial 
scales  to  provide  big  cities  with  fresh  water.  Small  villages  or 
settlements  in  rural  remote  areas  without  infrastructure  do  not 
profit  from  these  techniques.  The  technical  complexity  of  the  large 
plants  is  very  high  and  cannot  easily  be  scaled  down  to  very  small 
systems  and  water  demands.  Furthermore,  the  lack  of  energy 
sources  as  well  as  a  missing  connection  to  the  grid  complicates  the 
use  of  standard  desalination  techniques  in  these  places.  In  arid  and 
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Nomenclature 

Tic 

temperature  at  the  air-water  interface  in  the 

A 

air-water  exchanger  area  in  the  condensation  tower 

TP 

condensation  tower  (K) 
wall  temperature  (K) 

(m2) 

S 

water  solar  collector  area  (m2) 

a 

air-water  exchanger  area  in  the  evaporation  tower 

U 

overall  heat  transfer  coefficient  in  the  condensation 

Ce 

(m2) 

water  specific  heat  (J/kg/K) 

Uw 

tower  (W/m2/K) 

overall  energy  loss  from  the  absorber  to  outside 

Cw 

water  specific  heat  (J/kg/K) 

(W/m2/I<) 

Q 

moist  air  specific  heat  in  the  evaporation  tower 

w 

air  humidity  (kg  water/kg  dry  air) 

(J/kg/K) 

Wi 

saturation  humidity  in  the  evaporator 

cc 

moist  air  specific  heat  in  the  condensation  tower 

(kg  water/kg  dry  air) 

(J/kg/K) 

Wic 

saturation  humidity  in  the  condenser 

De 

water  mass  velocity  in  the  condensation  tower 

(kg  water/kg  dry  air) 

(kg/m2/s) 

z 

coordinate  in  the  flow  direction  (m) 

Dh  i 

hydraulic  diameter  of  air  flow  (m) 

X 

coordinate  in  the  flow  direction  (m) 

Dh2 

e 

hydraulic  diameter  of  water  flow  (m) 
thickness  of  the  condenser  plate  (m) 

Greek 

s 

c 

gravitational  acceleration  (m/s2) 
mass  velocity  of  humid  air  (kg/m2/s) 

a 

absorptance  of  the  collector  absorber  surface 

Gr 

Grashof  number 

Xq 

latent  heat  of  water  evaporation  (J/kg) 

hc 

heat  transfer  coefficient  for  the  condensate  film 

Ap 

wall  thermal  conductivity  (W/m/K) 

(W/m2/I<) 

Ae 

water  thermal  conductivity  (W/m/I<) 

hg 

air  heat  transfer  coefficient  at  the  air-water  interface 

Ac 

condensed  water  thermal  conductivity  (W/m/K) 

in  evaporator  (W/m2/I<) 

Ac 

humid  air  thermal  conductivity  in  the  condensation 

he 

heat  transfer  coefficient  at  the  water-condenser  inside 

tower  (W/m/K) 

wall  interface  (W/m2/I<) 

Pc 

water  density  (kg/m3) 

hi 

water  heat  transfer  coefficient  at  the  air-water 

Pc 

dynamic  viscosity  of  condensed  water  (Ns/m2) 

interface  in  evaporator  (W/m2/I<) 

T 

transmittance 

air  film  heat  transfer  coefficient  in  the  condensation 

n 

efficiency 

1 

tower  (W/m2/I<) 

width  of  water  solar  collector  (m) 

Subscripts 

L 

water  mass  flow  rate  in  the  evaporation  tower 
(kg/m2/s) 

1 

tower  bottom 

I 

solar  flux  (W/m2) 

2 

tower  top 

k 

water  vapor  mass  transfer  coefficient  at  the  air-water 

amb 

ambient 

interface  (W/m2/I<) 

e 

cooling  water  at  the  condensation  tower 

mw 

mc 

water  flow  rate  in  the  water  solar  collector  (kg/s) 
fresh  water  production  (kg/s) 

w 

water 

Pr 

Prandtl  number 

Re 

Reynolds  number 

g 

air  flowing  inside  the  evaporation  tower 

Pi 

T 

saturation  pressure  (Pa) 

G 

air  flowing  inside  the  condensation  tower 

temperature  (1<) 

l 

liquid 

Ti 

temperature  at  the  air-water  interface  in  the 
evaporation  tower  (K) 

semi  arid  regions  the  lack  of  drinkable  water  often  corresponds 
with  a  high  solar  insolation.  This  speaks  for  the  use  of  solar  energy 
as  the  driving  force  for  water  desalination  systems.  Especially  in 
remote  rural  areas  with  low  infrastructure  and  no  grid  connection, 
stand  alone  operating  systems  for  the  desalination  of  brackish 
or  sea  water  are  suitable  to  provide  small  settlements  with  clean 
potable  water. 

Solar  energy  has  been  receiving  greater  attention  in  recent  years 
for  various  applications  using  different  desalination  techniques.  Solar 
desalination  systems  can  be  divided  in  two  main  types,  i.e.  direct  and 
indirect  collection  systems.  The  direct  collection  systems  use  solar 
energy  to  produce  distillate  directly,  i.e.  in  the  solar  still,  which  has  a 
single  structure  and  does  not  require  sophisticated  technical  con¬ 
struction  or  operation  procedure.  The  solar  still  distillers  cannot 
compete  with  the  conventional  methods  such  as  multi-stage-flash 
(MSF),  multi-effect-distillation  (MED)  and  reverse  osmosis  (RO),  in 
particular  for  the  production  of  large  quantities  of  water.  However, 
they  seem  to  be  useful  in  providing  small  communities  with  fresh 


water  with  low  efficiency.  The  low  efficiency  is  originated  from  the 
fact  that  only  one  basin  is  used  in  which  the  brackish  water  is  heated, 
evaporated  and  then  condensed.  Many  solar  stills  have  been  con¬ 
structed  all  over  the  world.  To  mention  few  of  them,  we  cite  the 
works  of  Garg  and  Mann  [1],  Rajvanshi  [2],  Zaki  et  al.  [3],  Delyannis 
[4]  and  Baum  and  Bairamov  [5],  To  be  more  competitive  than  the 
classical  solar  desalination  system  as  solar  still,  others  promising 
solar  desalination  techniques  using  the  principle  of  humidification- 
dehumidification  (HD)  have  been  developed  for  the  production  of 
fresh  water  in  remote  and  sunny  regions  [6,7]. 

The  solar  powered  water  desalination  unit  based  on  the  HD 
principle  can  be  considered  as  one  of  the  performed  techniques 
developed  in  solar  desalination  field.  It  is  an  innovative  technology 
with  promising  diffusion  and  application  due  to  its  flexibility, 
simplified  design,  low  maintenance,  extended  life  time  for  over 
twenty  years,  quasi  null  energy  consumption,  low  capital  cost, 
construction  and  adaptation  for  use  in  rural  areas  to  produce  fresh 
water  for  drinking  and  irrigation  [8-10], 
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As  mentioned  previously,  the  HD  process  is  an  interesting  techni¬ 
que  that  has  been  widely  adapted  for  water  desalination.  However, 
the  major  inconveniency  of  solar  desalination  unit  especially  the  HD 
technique  is  their  limited  fresh  water  production.  Therefore,  the 
improvement  of  small  solar  desalination  units  presents  an  attractive 
solution  to  decentralised  water  demand,  which  is  the  case  of  remote 
and  rural  area  in  developing  countries.  There  are  a  number  of  studies 
on  the  humidification-dehumidification  desalination  process  most  of 
which  focus  on  performance  evaluation  and  efficiency  improvement 
in  order  to  optimize  the  fresh  water  production  of  solar  desalination 
units  and  to  reduce  their  costs. 

Veera  et  al.  [11]  studied  the  feasibility  of  round  the  clock  operation 
of  a  low  temperature  desalination  system  powered  by  solar  collectors 
with  the  use  of  thermal  energy  storage.  Results  from  the  theoretical 
study  confirm  that  thermal  energy  storage  is  a  useful  component  of 
the  system  for  conserving  thermal  energy  to  meet  the  energy  demand 
when  direct  solar  energy  resource  is  not  available.  It  has  been 
determined  that  a  solar  collector  area  of  18  m2  with  a  thermal  energy 
storage  volume  of  3  m3  is  adequate  to  produce  100  L/d  of  freshwater 
round  the  clock  considering  fluctuations  in  the  weather  conditions. 

Karan  et  al.  [12]  used  multi-parameter  optimization  methods  to 
find  the  optimal  operating  conditions  for  both  CAOW  (Closed  air  open 
water)  and  OAOW  (open  air  open  water)  humidification-dehumidifi- 
cation  desalination  cycles  that  result  in  maximum  gained  output  ratio 
(GOR)  and  to  investigate  some  characteristics  of  the  various  cycles,  as 
well  as  open  air  open  water  cycles,  each  with  either  an  air  or  a  water 
heater,  were  considered  in  this  analysis.  Numerical  optimization 
resulted  in  a  substantial  increase  in  GOR  for  all  four  cycle  types 
compared  to  previous  best-case  conditions  found  using  heuristic 
studies.  The  GOR  of  the  cycles  was  found  to  decrease  with  increasing 
component  terminal  temperature  difference  (TTD).  In  addition,  differ¬ 
ent  cycles  perform  best  at  different  temperature  differences. 

In  the  same  year,  Mehrgoo  and  Amidpour  [13]  developed  a 
model  and  a  structured  procedure  to  optimize  the  shape  and 
structure  (relative  sizes,  aspect  ratios)  of  a  multi-effect  humidifi- 
cation-dehumidification  (MEH)  desalination  unit.  The  results 
showed  that  the  inlet  cold  and  hot  water  temperatures  and  the 
column  heights  play  important  roles  in  the  constructal  design  of  a 


MEH  unit,  especially  as  the  total  volume  increases.  Moreover,  the 
tube  diameters  and  the  number  of  condenser  tubes  should  be 
suitably  selected  based  on  the  values  of  the  condenser,  the 
humidifier  inlet  water  temperatures  and  the  column  heights. 

Narayan  et  al.  [14]  and  Mistry  et  al.  [15]  investigated  the 
thermodynamic  behavior  of  HD  cycles  using  the  First  and  Second 
Laws  of  Thermodynamics,  respectively.  Both  papers  discussed 
methods  for  improvement  and  optimization  of  the  cycles  based 
on  thermodynamic  arguments  and  single  parameter  optimization. 

Zamen  et  al.  [16]  optimized  the  total  solar  HD  system.  Such 
optimization  was  mainly  intended  to  reduce  the  fresh  water  produc¬ 
tion  costs.  They  showed  that  solution  obtained  by  cost  objective 
function  results  in  a  cost  7-28%  lower  than  other  objective  functions. 

The  present  work  deals  with  modeling,  simulation  and  experi¬ 
mental  validation  of  a  new  generation  of  water  desalination  unit 
by  solar  energy  using  a  multiple  effect  humidification  process.  The 
main  objectives  of  the  present  paper  are: 

•  To  present  a  steady  state  mathematical  model  of  the  different 
components  of  the  unit  (water  solar  collector,  evaporation 
tower  and  condensation  tower). 

•  To  investigate  the  thermal  performances  of  the  solar  desalina¬ 
tion  unit  exposed  to  a  variation  of  operating  key  parameters. 

•  To  present  numerical  and  experimental  data  as  an  example  of 
the  validation  process  that  has  been  carried  out  in  order  to 
assess  the  credibility  of  the  numerical  models. 

The  present  paper  is  organized  as  follows.  Sections  2  and  3  are 
dedicated,  respectively,  to  the  process  description  and  the  devel¬ 
oped  mathematical  models.  Section  4  presents  the  computa 
tion  process.  The  instrumentation  and  the  simulation  results  are 
presented,  respectively  in  Sections  5  and  6. 

2.  Process  description 

Figs.  1  and  2  show,  respectively,  a  schematic  diagram  and  a  photo 
of  the  experimental  setup  of  the  humidification  and  dehumidification 
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Fig.  1.  Schematic  diagram  of  the  HD  water  desalination  unit. 
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Fig.  2.  A  photo  of  the  solar  desalination  unit  experimental  setup. 
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Fig.  3.  Thermal  energy  balance  of  differential  section  of  solar  water  collector. 


(HD)  water  desalination  unit.  The  developed  unit  operates  at  atmo¬ 
spheric  pressure  using  air  as  a  carrier  for  vapour.  Three  compartments 
constitute  the  installation:  solar  collectors,  evaporation  tower  and 
condensation  tower. 

The  principle  of  functioning  of  this  desalination  process  is  as 
follows.  First,  brackish  water  or  sea  water  is  heated  by  the  solar 
collectors  as  much  as  possible.  Then,  the  hot  water  is  injected  to 
the  top  of  the  evaporation  tower.  A  pulverizer  with  a  special  shape 
is  used  to  assure  a  uniform  pulverisation  of  the  hot  water  in  all  the 
sections  of  the  tower.  Hot  and  saturated  air  mixes  with  the  rising 
air  current  toward  the  condensation  tower.  Then  it  condenses  in 
contact  with  the  cold  condensation  plates.  The  circulation  of  the 
air  in  the  evaporation  tower  may  occur  in  natural  convection  or  in 
forced  convection  with  two  functioning  modes:  closed  air  circuit 
or  open  air  circuit.  In  the  course  of  working  in  free  convection,  the 
displacement  of  humid  air  in  the  ascending  direction  of  the  tower 
is  caused  by  the  difference  of  the  density  and  temperature  of 
the  humid  air  in  the  tower.  The  forced  convection  is  insured  by  a 
helical  fan  fixed  at  the  bottom  or  at  the  top  of  the  tower.  The 
choice  of  the  air  circulation  mode  in  the  tower  depends  on  the 
shape  and  the  dimension  of  the  packed  bed. 

The  detailed  descriptions  of  the  unit  main  components  are  as 
follows: 

The  solar  collector  is  the  major  component  of  the  unit.  In  fact, 
its  characteristics  have  an  important  influence  on  the  operation 
and  efficiency  of  the  unit.  It  has  an  absorber  plate  containing 
uniformly  spaced  parallel  copper  pipes  with  its  upper  surface 


painted  matt  black  to  increase  the  absorptivity  of  the  system.  The 
absorber  can  be  also  made  from  a  polypropylene  material  with 
very  thin  and  tightly  spaced  capillary  tubes.  The  fluid  circulates  in 
a  forced  convection  and  in  one  direction  in  the  absorber,  which  is 
covered  with  a  single  glass  of  high  transmissivity  to  solar  radia¬ 
tion.  All  the  parts  are  fitted  inside  an  external  case.  The  undersides 
of  the  absorber  and  the  side  casing  are  well  insulated  to  reduce  the 
heat  losses  to  ambient  air. 

The  evaporation  tower  produces  the  water  vapour.  It  is 
equipped  with  packed  bed  to  increase  the  contact  surface  and 
therefore  improve  the  humidification  rate.  Due  to  the  nature  of 
the  water  (hot  salinity  level,  chalky,  solid  residuals  and  corrosion 
problems,  etc.),  the  packed  bed  used  in  the  unit  must  be  carefully 
chosen.  Thorn  trees  or  palm  tree  leaves  are  well  suited  for  this 
application.  They  are  abundant,  free  and  have  a  high  resistance 
against  the  forth-mentioned  problems. 

The  condensation  chamber  contains  polypropylene  condensa¬ 
tion  plates  through  which  the  cold  salty  water  circulates  for 
preheating.  The  structure  of  the  evaporation  tower  and  the 
condensation  tower  is  made  of  aluminium.  The  thermal  insulation 
is  achieved  by  Styrofoam  layers  with  a  thickness  of  40  to  50  mm. 
To  insure  isolation  against  vapour  and  water  circulating  in  the  unit 
and  to  protect  the  Styrofoam  plates  against  corrosion  caused  by 
brackish  water,  the  inside  of  the  unit  is  covered  with  polypropy¬ 
lene  plates.  The  distilled  water  is  collected  in  a  basin  at  the  bottom 
of  the  condensation  tower.  Water  not  evaporated  is  collected  in  a 
basin  at  the  bottom  of  the  evaporation  tower  and  then  recycled  or 
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discarded  in  the  case  that  the  salt  rate  is  too  high.  The  two  basins 
are  made  of  polypropylene  material.  Pumps  are  installed  to 
facilitate  water  circulation  in  the  different  parts  of  the  desalination 
unit.  Polypropylene  ducts  are  used  due  to  their  resistance  to 
corrosion. 


3.  Process  modeling 


3.2.  Evaporation  tower  mathematical  model 

The  model  of  the  evaporation  tower  is  a  set  of  equations 
developed  using  thermal  and  mass  balances  for  the  water  phase, 
the  air  phase,  and  the  air-water  interface  on  an  element  of  volume 
of  height  dz  (Fig.  4).  The  approach  is  similar  to  that  used  in  the 
cooling  tower  study  [17,18]. 

To  establish  the  model  for  the  system,  the  following  assump¬ 
tions  are  made: 


Mathematical  modeling  is  at  the  very  core  of  the  design  and 
optimization  process  for  thermal  systems  because  the  mathema¬ 
tical  model  brings  out  the  important  variables  and  the  governing 
parameters.  We  have  developed  a  mathematical  model  in  steady 
state  regime  based  on  heat  and  mass  transfers  in  each  component 
of  the  unit  i.e.  water  solar  collector,  evaporation  tower  and 
condensation  tower  of  the  desalination  unit  in  order  to  numeri¬ 
cally  simulate  the  HD  system.  The  developed  models  allow  a  good 
description  of  the  real  process  and  a  better  monitoring  of  the 
spatial  distribution  of  the  different  parameters  like  air  and  water 
temperatures  and  air  humidity. 


3.1.  Water  solar  collector  mathematical  model 

Since  there  are  no  mass  or  concentration  changes,  the  water 
solar  collector  can  be  modeled  solely  from  thermal  energy  bal¬ 
ances.  To  establish  the  model  for  the  system,  the  following 
assumptions  are  considered: 

-  The  speed  of  the  fluid  is  uniform,  therefore  the  local  state  of  the 
fluid  and  of  the  absorber  depends  only  on  one  side  x  that  varies 
from  0  to  /. 

-  The  absorber  and  fluid  have  the  same  temperature  at  any  point. 

-  The  fluid  temperature  remains  under  the  100  °C  point. 

The  thermal  energy  balance  equations  for  the  system  formed 
by  the  absorber  and  water  for  a  slice  of  the  collector  with  a  width 
of  /,  a  length  of  dx  and  a  surface  of  ds  as  shown  in  Fig.  3,  is  the 
next  one: 

[(The  quantity  of  energy  received  by  a  surface  absorber  ele¬ 
ment,  ds]  =  [(the  quantity  of  energy  gained  by  water)+(the  overall 
losses  toward  the  outside)] 

ha  ds  =  Qyv  +  Q\ossab_amb  ( 1 ) 

with 

Qw  =  mwCwdTw 


Q_lossab_amb  —  UW(TW  Tamb)  ds 

Eq.(l)  can  be  written  as: 

I  t  a  ds  =  mwCwdTw  +  UW(TW—Tamb)ds 

It  can  be  also  further  simplified  as: 

dTw  Uwl  (zal  |  T  \ 
dx  ~  mwCw  \UW  amb  WJ 


O') 


(1") 


The  mathematical  model  for  the  solar  collector,  allows  the 
study  for  constant  insolation  level  and  constant  ambient  tempera¬ 
ture,  the  fluid  temperature  variation  at  the  outlet  of  the  collector 
as  a  function  of  the  control,  geometrical,  and  physical  parameters 
(flow  rate,  collector  area,  material,  and  inclination  angle  of  the 
collector).  In  addition,  the  fluid  temperature  and  flow  rate  at  the 
collector  outlet  are  the  two  parameters  with  the  most  significant 
impact  on  the  unit  production  as  they  are  among  the  input 
parameters  of  the  evaporation  tower. 


-  The  tower  is  adiabatic. 

-  Head  losses  in  the  tower  column  are  low  and  it  can  be  assumed 
that  the  density  of  the  air  mass  flux  is  constant. 

-  Since  the  evaporation  rate  is  small  compared  to  the  water  flow 
rate,  then  the  water  mass  flux  remains  constant  throughout  the 
column. 

-  The  water  distribution  over  the  packed  bed  is  uniform. 

•  For  the  water  phase 

The  water  thermal  balance  in  an  element  of  volume  of  the 
evaporation  tower  of  height  dz,  is  given  by  the  following 
equation: 

[The  quantity  of  heat  carried  by  water]  =  [the  quantity  of  heat 
transmitted  toward  the  surface  of  separation  water-air] 
Mathematically  this  can  be  expressed  as: 

LCidTi  =  hia(Ti~Tj)dz  (2) 


It  can  be  further  simplified  as: 
dT,  h,a(T,-Ti) 


dz 


LC, 


(2) 


L 

Tl2 


L 

Ti 
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Fig.  4.  An  element  of  volume  of  the  evaporation  chamber  of  height  dz. 
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•  For  the  air  phase 

The  air  thermal  balance  in  an  element  of  volume  of  the 
evaporation  tower  of  height  dz,  is  the  following  one: 

[The  quantity  of  heat  received  by  the  surface  of  separation 
water-air]  =  [the  quantity  of  heat  transmitted  toward  the  air 
flow  rate] 

Mathematically  this  can  be  expressed  as: 


beforehand.  To  do  so,  a  set  of  experiments  were  carried  out  on  a 
test  bench  for  a  cooling  tower  to  determine  the  different  specific 
heat  and  mass  exchange  coefficients  at  the  contact  surface 
between  the  garnishing  and  the  humid  air  [9[.  The  relationships 
giving  the  mass  exchange  coefficients  ( kg )  and  the  heat  exchange 
coefficients  hi  as  a  function  of  the  air  flow  rate  (G)  and  water  flow 
rate  (L)  are 


hga(Tj-Tg)dz  =  GCgdTg 


(3) 


kg  = 


2.09  G 


■0.11515,0.45 


(8) 


It  can  be  further  simplified  as: 

dTg  _  hga(Tj-Tg) 
dz  ~  GCg 


(3') 


•  Air-water  interface 

At  the  air-water  interface,  the  thermal  balance  can  be  written 
under  the  following  shape: 

[The  quantity  of  heat  transmitted  from  the  water  current  through 
the  liquid  film  toward  the  surface  of  separation  water-air]  =  [ 
(the  quantity  of  heat  transmitted  from  the  surface  of  separation 
through  the  air  film  toward  the  air  current)+(the  heat  quantity 
provided  to  evaporate  the  water  mass  quantity  transferred  of  the 
liquid,  through  the  interface  toward  the  air  current)] 


Mathematically  this  can  be  expressed  as: 
hia(T,-Ti)  =  hga(Tj-Tg)  +  20kga(W  ,-Wg) 

It  can  be  further  simplified  as: 
w.-w  hia(T,-Ti)  +  hga(Tg-Ti) 


(4) 

(4') 


Also,  it  can  be  expressed  as  [19]: 

Wi  =  0.62198  (5) 

where  Pws  is  the  saturation  pressure  of  the  water  vapor  corre¬ 
sponding  to  a  certain  temperature  is  computed  from  [20]  as: 

InfPws)  =  -6096.938  =1  +  21 .240964-2.71 1 191 0“2T, 

'  i 

+  1.6739510“5Tf  +  2.433501n(T;)  (6) 


At  the  air-water  interface,  the  mass  balance  is  given  by  the 
following  equation: 

[The  mass  water  quantity  transferred  from  the  interface  through 
the  air  film  toward  air  flow  rate]  =  [the  quantity  of  steam  transmitted 
to  air  flow  rate] 

Mathematically  this  can  be  expressed  as: 
kga(W~Wg)dz  =  GdWg  (7) 

It  can  be  further  simplified  as: 
dWg  _  kgajWt-Wg) 

dz  ~  G  1 

The  steady  state  model  of  the  evaporation  tower  using  the 
equations  presented  above  could  be  used  to  monitor  the  influence 
of  the  command  parameters  on  the  tower  operation.  In  particular, 
it  allows  to  determine  the  temperature  and  the  amount  of  water  in 
the  air  at  the  outlet  of  the  evaporation  tower  as  a  function  of  the 
water  temperature  at  the  entrance  (top  of  the  tower)  for  a  given 
tower  height.  The  obtained  model  can  be  used  to  size  this  type  of 
evaporation  tower  for  desired  water  and  air  temperatures. 

Due  to  the  nature  of  water  (hard  salinity,  deposits  of  lime 
stones,  solid  residues,  corrosion  problems,  etc.),  the  garnishing  to 
be  used  in  this  application  must  be  chosen  carefully.  For  instance, 
the  thorn  trees,  abundant  and  free  are  appropriate  under  these 
operating  conditions  and  their  characteristics  need  to  be  specified 


,  5900  G 

h,= - 


0.5894^0.169 

a 


(9) 


The  air-film  heat  transfer  coefficient  and  the  mass  transfer 
coefficient  on  the  air-water  interface  are  coupled  by  the  Lewis 
relation  [21]  as  follows: 

hg=  Ggkg  (10) 


3.3.  Condensation  tower  mathematical  model 

The  formulation  of  the  mathematical  model  of  the  condensa¬ 
tion  tower  in  the  steady  state  regime  is  based  on  thermal  and  mass 
balances.  It  allows  developing  the  coupling  equations  between  the 
temperature  of  the  cooling  water  of  the  condenser  and  the  humid 
air  temperature  and  water  content.  The  condensation  rate  is 
determined  using  an  algebraic  equation  that  relates  the  variation 
of  the  water  content  with  the  height  of  the  tower.  The  balances  are 
done  on  an  element  of  volume  of  the  tower  of  height  dz  (Fig.  5). 

To  establish  the  model  for  the  system,  the  following  assump¬ 
tions  are  made: 

-  The  tower  is  adiabatic. 

-  Each  stage  is  completely  homogeneous  as  to  avoid  thermal 
gradients  in  the  direction  normal  to  the  condenser  plates. 

-  The  length  dz  of  the  tube  element  is  small;  therefore,  the 
temperature  in  each  stage  changes  linearly  throughout  the 
length  of  the  element. 

-  Colburn  and  Flougen  method  developed  by  Votta,  which  states 
that  the  sensitive  heat  transfer  of  the  condensed  water  droplets 
is  negligible,  is  used  [22]. 

-  The  air  heat  and  mass  transfer  coefficients  at  the  air-water 
interface  are  related  by  the  Lewis  equation. 

-  The  condensate  film  flow  is  normally  streamlined. 

-  The  heat  flows  through  the  film  by  conduction. 
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Fig.  5.  An  element  of  volume  of  the  condensation  chamber  of  height  dz. 
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The  mathematical  model  for  the  condensation  tower  is  formed 
by  the  following  system  of  equations: 

•  For  the  water  phase 

[The  quantity  of  heat  transmitted  from  air-water  interface  to 
water  in  the  condenser]  =  [the  quantity  of  heat  carried  away 
by  water] 

Mathematically  this  can  be  expressed  as: 

UA(Tic-Te)dz  =  DeCedTe 

It  can  be  further  simplified  as: 
dTe  UA(Tic-Te) 
dz  DeCe 

•  For  the  air  phase 
[The  quantity  of  heat  carried  by  air]  =  [(sensible  heat  quantity 
transmitted  toward  the  surface  of  separation  air-water)-(latent 
heat  quantity  due  to  the  condensation  of  steam  transmitted  to  the 
surface  of  separation  air-water)] 

Mathematically  this  can  be  expressed  as: 

GCGdTG  =  hGA(TG-TjC)dz—A0kGA(WgC-WjC)dz  (12) 

It  can  be  further  simplified  as: 

dTp  hG(TG-TjC )  +  A0kG(WG-WjC))  (12') 

dz  ~  CCG 


(ID 


(11') 


For  a  flow  between  two  vertical  plates,  the  heat  transfer 
coefficient  hc  is  given  by  [17]: 

hc=  ^-0.779(GrG4)  (18) 

Vh\ 

The  coefficient  kG  is  determined  using  the  Lewis  relation: 

kc=p-  (19) 

<-G 


The  overall  heat  transfer  coefficient  from  the  air-condensate 
interface  to  the  cooling  water  inside  the  condenser  is  approxi¬ 
mated  by  [9[: 


where  he  is  given  by  [23]: 

he  =  0.023  Re08Pr°  33 
Uh2 


and  hc  is  expressed  using  the  Nusselt  relation  [23]: 


,4/  P2cS^3c 

V  4/Jcz(TjC—Tp) 


(21) 


(22) 


The  steady  state  model  developed  allows  finding  the  influence 
of  the  command  parameters  on  the  tower  operation  and  the  sizing 
of  this  type  of  condensation  tower  for  desired  amount  of  distilled 
water  mc  produced  by  the  unit. 


•  At  the  air-condensate  interface 
The  heat  balance  at  air-condensate  interface  is  given  by  the 
following  equation: 

[(cooling  sensible  heat  quantity  of  humid  air  received  by  air- 
condensate  interface )+( latent  heat  quantity  liberated  during 
the  condensation  at  air-condensate  interface)]  =  [quantity  of  heat 
transmitted  from  air-condensate  interface  to  the  cooling  water] 
Mathematically  this  can  be  expressed  as: 

hGA(TG-Tic)dz  +  X0kcA(Wc-W  ic)dz  =  UA(Tic-Te)dz  (13) 


It  can  be  further  simplified  as: 
Wic  =  Wc  +  ^cCTc-fic)  +  U(Te-Tic) 

20kc 

also,  it  can  be  expressed  as  [19]: 


(13') 


4.  Computational  process 

The  model  equations  developed  in  the  preceding  section  are 
composed  of  a  set  of  highly  interlinked  nonlinear  equations. 
Therefore,  iterative  solution  is  necessary  to  calculate  the  system 
characteristics  for  each  unit  components.  Moreover,  since  these 
components  are  coupled  and  the  output  of  each  component  is  the 
input  of  the  next,  it  is  necessaiy  to  couple  the  various  components 
of  the  present  desalination  unit  and  get  a  global  steady  state 
mathematical  model  that  would  provide  a  good  description  of  the 
entire  system.  Fig.  6  gives  a  general  picture  of  the  flow  chart  of  the 
computer  model. 

A  computation  program  was  compiled  in  terms  of  the  system 
model  to  predict  the  system  performance.  The  calculation  proce¬ 
dure  of  the  program  is  as  follows: 


Wic  =  0.62198^  (14) 

where  P,  is  the  saturation  pressure. 

The  mass  balance  at  air-condensate  interface  is  given  by  the 
following  equation: 

[The  quantity  of  air  transferred  toward  air-water  interface]  = 
[the  quantity  of  steam  condensed] 

Mathematically  this  can  be  expressed  as: 

I<  CA  ( Wgc  -WjC)dz  =  CdWG  (15) 


It  can  be  further  simplified  as: 


dWG 

dz 


=  kGA 


( WG-Wic ) 

c 


(15') 


Following  the  water  balance  equation: 
dmc  =  CdWG 


(16) 


So  the  flow  rate  of  the  condensed  water  is: 


mw,  T„ 


- Water  loop 


dmc  =  kcA(W  ic-WG)dz 


(17) 


Fig.  6.  A  general  picture  of  the  flow  chart  of  the  computer  model. 
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1.  Input  the  known  data  for  each  unit  components,  such  as  the 
functioning  parameters  and  the  geometric  parameters. 

2.  Calculate  the  data  concerning  the  water  solar  collector  and 
obtain  the  temperature  of  feed  water  to  the  evaporator. 

3.  Calculate  the  data  concerning  the  evaporation  tower  and  obtain 
evaporation  rate,  outlet  temperature  of  seawater,  outlet  tem¬ 
perature  and  humidity  of  the  air. 

4.  Calculate  the  data  concerning  condenser  to  obtain  the  outlet 
temperature  and  humidity  of  the  air,  outlet  temperature  of  the 
cooling  water  and  the  fresh  water  production  in  particular. 

5.  Calculate  the  indexes  of  system  performance  and  output  the 
results. 


5.  Instrumentation 

The  measured  variables  in  the  experiment  include  inlet  and 
outlet  air  and  water  temperatures  in  each  unit  components, 
ambient  temperature,  air  velocity  and  the  solar  irradiation  inci¬ 
dent  on  the  field  collector's  plane.  The  water  solar  collector  was 
instrumented  with  Pt  100  thermistors  for  measuring  the  outlet 
and  the  inlet  of  water  temperatures  and  ambient  temperature.  The 
Pt  100,  which  measured  the  ambient  temperature,  was  kept  in  a 
shelter  to  protect  the  sensor  from  direct  sunlight.  The  inlet  and  the 
outlet  water  and  air  temperatures  in  both  the  evaporation  tower 
and  the  condensation  tower  were  measured  using  the  TH  100 
transmitter.  The  water  inlet  flow  rate  was  measured  using  a  flow 
meter,  while  the  fresh  water  flow  rate  was  measured  by  graduated 
cylinder.  A  pyranometer  was  used  to  measure  the  total  solar 
irradiation  placed  in  a  horizontal  plane  adjacent  to  the  collector. 
All  sensors,  which  were  calibrated  before  using  to  determine  the 
probes  sensibility,  were  connected  to  a  data  acquisition  system 
(type  Agilent  34970  A).  The  detailed  technical  specifications  of 
sensors  and  probes  used  in  experimental  set  up  are  presented  in 
Table  1. 


the  solar  collector  by  internally  recovering  the  latent  heat  of 
vaporisation.  For  this  purpose,  the  feed  water  was  preaheated  in 
the  condenser  by  exchange  of  heat  with  the  vapour  at  condensa¬ 
tion.  The  coolant  water  passes  the  condenser  plates  in  counter 
flow  to  the  hot  feed  water  and  is  gradually  warmed  by  the  latent 
heat  of  condensation.  When  the  coolant  water  comes  out  of 
the  condenser  it  is  fed  to  the  flat  plate  solar  collector  for  further 
heating.  In  thermal  desalination  processes  the  gained  output  ratio 
(GOR)  and  the  thermal  efficiency  of  the  solar  collector  are  the  two 
most  commonly  encountered  parameters  for  the  evaluation  and 
assessment  of  such  systems.  The  GOR  can  either  be  defined  as  a 
mass  ratio  (kg  distillate  produced/kg  steam  consumed)  or  as  an 
energy  ratio.  The  drawback  of  the  mass  ratio  definition  is  that  it 
ignores  the  amount  of  heat  transferred  from  the  steam  to  the 
product  water  and  reject  brine.  As  an  energy  ratio  the  GOR  is 
defined  as  the  total  latent  heat  of  evaporation  of  the  product  water 
to  the  input  thermal  energy.  Theoretical  energy  required  to 
produce  the  distillate  divided  by  the  actual  thermal  energy 
consumed  in  the  evaporator.  Mathematically,  the  GOR  of  the  solar 
desalination  unit  can  be  written  as  [24,25]: 


GOR  = 


mcA0 

LCi(Ta-Tn) 


(23) 


where  mc  is  the  condensation  flow  rate,  L  is  the  water  mass 
flow  rate  in  the  evaporation  tower,  and  Aa  is  the  latent  heat  of 
vaporization. 

The  thermal  efficiency  of  the  water  solar  collector  was  eval¬ 
uated  by  comparing  the  total  radiant  heat  energy  flux  to  the 
energy  gained  water  after  passage  through  the  collector.  This  is 
given  in  the  following  equation  [24]: 

mwCw(Two-Twi) 

"= - Is -  (24) 

where  mw  is  the  water  mass  flow  rate,  Cw  is  the  water  heat 
capacity,  Two  is  the  water  outlet  temperature,  Tw,  is  the  water  inlet 


6.  Results  and  discussion 


Table  2 

Characteristics  of  the  pilot  desalination  unit. 


Simulation  is  the  process  of  subjecting  the  model  for  a  given 
thermal  system  to  various  inputs,  such  as  operating  conditions,  to 
determine  how  it  behaves  and  thus  predict  the  characteristics  of 
the  actual  physical  system.  Though  simulation  may  be  carried  out 
with  scale  models  and  prototypes,  the  expense  and  effort  involved 
generally  makes  it  impossible  to  use  these  for  design  because 
many  different  designs  and  operating  conditions  need  to  be 
considered  and  evaluated.The  characteristics  of  the  pilot  desalina¬ 
tion  unit  used  in  simulation  are  listed  in  Table  2. 

6.1.  Thermal  performance 

The  solar  desalination  unit  was  designed  originally  to  reduce 
the  thermal  energy  consumption  necessary  for  heating  water  in 


Solar  collector 

Area 

Effective  transmission  absorption 
Riser  tube  material 
Absorber  surface 
Loss  coefficient 

Back  insulation,  thickness,  mm 

Evaporation  chamber 

Size 

Packed  bed  type 
Mass  transfer  coefficient 
Heat  transfer  coefficient 

Condensation  tower 

Size 

Plates  type 


7.20  m2 
0.83 
Copper 

Pain  mat  black 
3.73  W/m2/I< 

Fibre  glass,  50 

1.20  m  x  0.50  m  x  2.55  m 
Solid  packing:  thorn  trees 
kg=( 2.09  x  G011515  x  L045)/a 
hi=(5900  x  G05894  x  L0A69)/ahg=Cg  x  kg 

1.2  m  x  0.36  m  x  3.0  m 
Polypropylene 


Table  1 

Sensors  and  probes  used  in  experimental  set  up. 


Sensor/type 

Reference 

Description 

Value 

Pt  100 

PRO-SL  100  PV 

Sensibility 

0.3799  £2/°C 

TH  100 

KIMO  TH  100-AOD 

Humidity 

Accuracy 

±2% 

Sensibility 

0.159  mA/% 

Temperature 

Accuracy 

±  0.2  C 

Sensibility 

0.160  mA/°C 

Pyranometer  CE  180 

CIMEL  CE  180 

Accuracy 

1% 

Sensibility 

12.29  nV/Wm2 
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temperature,  I  is  the  solar  irradiation  intensity  and  S  is  the 
collector  area. 

Figs.  7  and  8  show  the  thermal  efficiency  variation  of  water 
solar  collector  as  a  function,  respectively,  of  water  mass  flow  rate 
and  inlet  water  temperature  for  different  values  of  solar  radiation. 
Fig.  7  presents  that  the  thermal  efficiency  of  the  water  solar 
collector  increases  with  the  increase  of  the  water  flow  rate.  This 
may  be  due  to  the  fact  that  for  high  flow  rates,  the  collector 
operating  temperature  would  be  lower,  resulting  in  lower  heat 
losses  and,  subsequently,  higher  efficiency.  In  addition,  for  water 
flow  rate  values  higher  than  0.4  kg/s  one  notes  that  an  asymptotic 
regime  was  reached  during  which  the  increase  in  the  water  flow 
rate  is  useless.  Fig.  8  shows  that  the  thermal  efficiency  of  the  water 
solar  collector  decreases  with  an  increase  of  the  inlet  water 
temperature  (Twi).  On  one  hand,  for  inlet  water  temperature  values 
inferior  to  ambient  temperature,  the  thermal  efficiency  increases 
with  a  decrease  of  solar  radiation  intensity.  On  the  other  hand,  for 
inlet  water  temperature  values  superior  to  ambient  temperature, 
the  thermal  efficiency  decreases  with  a  decrease  of  solar  radiation 
intensity.  From  Fig.  8,  one  can  also  notes  that  for  an  inlet  water 
temperature  equal  to  the  ambient  temperature,  the  solar  radiation 
intensity  has  no  effect  on  the  thermal  efficiency  of  the  water  solar 
collector. 

The  variation  of  the  gain  output  ratio  with  cooling  water  mass 
flow  rate  at  the  condensation  tower  and  water  mass  flow  rate  at 
the  evaporation  tower  are  presented,  respectively  in  Figs.  9  and  10. 
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Fig.  7.  Effect  of  the  water  mass  flow  rate  on  the  collector  thermal  efficiency. 
(Twi=20  “C,  Tamb=30  °C). 


Fig.  8.  Effect  of  the  inlet  water  temperature  on  the  collector  thermal  effeciency  for 
different  values  of  solar  radiation  ( mw — 0.1  kg/s,  Tamb= 30  C). 


A  slight  decrease  of  the  GOR  is  observed  in  Fig.  9.  When  cooling 
water  mass  flow  rate  is  increased  and  then  stabilizes.  Fig.  10 
illustrates  that  the  optimal  water  mass  flow  rate  is  about  0.4  kg/s, 
of  which  the  GOR  is  maximum.  The  GOR  of  the  system  is  3.  This 
result  can  be  explained  by  the  fact  that  while  increasing  the  water 
mass  flow  rate  inside  the  evaporation  tower,  the  rate  of  vaporisation 
increases  which  in  its  part  augments  the  fresh  water  production, 
hence,  the  GOR. 

6.2.  Experimental  validation 

The  validation  of  the  model  developed  for  a  given  system 
is  another  very  important  consideration  because  it  determines 
whether  the  model  is  a  faithful  representation  of  the  actual 
physical  system  and  indicates  the  level  of  accuracy  that  may  be 
expected  in  the  predictions  obtained  from  the  model.  To  validate 
the  developed  models  of  the  desalination  unit,  a  series  of  experi¬ 
ments  was  conducted  using  the  pilot  desalination  unit  located  at 
the  National  Engineering  School  of  Sfax  (E.N.l.S)— Sfax  University, 
Tunisia.  Experimental  measurements  were  taken  on  the  solar 
collector  and  the  distillation  module  (evaporation  and  condensa¬ 
tion  chambers). 

Fig.  11  illustrates  a  comparison  between  predicted  and  mea¬ 
sured  collector  outlet  temperature  (Two)  as  a  function  of  solar 
radiation.  It  is  apparent,  as  expected,  that  for  a  constant  water  flow 


Fig.  9.  Effect  of  the  cooling  water  mass  flow  rate  at  the  condensation  tower  on  GOR 
(Tei  =  20  °C,  L=0.2  kg/s,  G=0.1  kg/s). 


Fig.  10.  Effect  of  the  water  mass  flow  rate  at  the  evaporation  tower  on  GOR 
(Tei  =  20  °C,  De=0.1  kg/s,  G=0.1  kg/s). 
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rate  the  outlet  water  temperature  increases  with  solar  radiation.  It 
can  be  also  seen  that  a  good  agreement  had  been  found  between 
the  simulated  and  the  experimental  results. 

Fig.  12  shows  a  comparison  between  the  experimental  and  the 
theoretical  values  of  air  absolute  humidity  ( Wg2 )  as  a  function  of 
inlet  water  temperature  at  the  level  of  the  evaporation  tower.  It  is 
clear  from  this  figure  that  the  amount  of  water  in  the  humid  air  at 
the  top  of  the  evaporation  tower  increases  significantly  as  the 
water  temperature  (T(2)  at  the  evaporation  tower  inlet  increases.  It 
can  be  also  noticed  that  the  theoretical  curve  of  water  content  Wg2 
as  a  function  of  the  temperature  Ti2  exceeds  the  experimental 
curve,  as  seen  in  Fig.  12.  The  gap  between  the  theoretical  and 
experimental  values  increases  with  T(2.  This  may  be  due  to  the 
thermal  losses  that  increase  with  the  temperature,  even  though 
the  process  is  assumed  to  be  theoretically  adiabatic. 

Fig.  13  shows  a  comparison  between  the  measured  and  the 
calculated  cooling  water  temperature  at  the  exit  of  the  condensa¬ 
tion  tower  as  a  function  of  inlet  water  temperature  at  the  entrance 
of  the  evaporation  tower.  Through  this  figure  one  can  note  that 
increasing  the  evaporation  tower  inlet  water  tempertaure  leads  to 
a  significant  increase  in  the  condensation  tower  outlet  cooling 
water  temperature.  This  may  be  due  to  the  fact  that  an  increase  of 
the  inlet  water  temperature  in  the  evporator  increases  the  outlet 
moist  air  temperature  leaving  the  latter  to  the  condenser.  Conse¬ 
quently,  the  outlet  cooling  water  temperature  increases  due  to  the 


latent  heat  of  condensation.  It  was  also  noticed  that  the  experi¬ 
mental  values  were  less  than  the  theoretical  ones.  This  may  be 
attributed  partly  to  the  fact  that  the  theoretical  model  assumes 
that  each  stage  is  completely  homogeneous  as  to  avoid  thermal 
gradients  in  the  direction  normal  to  the  condenser  plates. 

Figs.  14  and  15  show  a  comparison  between  the  experimental 
and  the  theoretical  values  of  the  outlet  fresh  water  flow  rate  (mc) 
as  a  function,  respectively,  of  the  water  temperature  at  the 
evaporator  inlet  and  the  air  temperature  at  the  condenser  inlet. 
From  this  figure  we  can  note  that  the  flow  rate  of  the  fresh  water 
increases  with  increasing  air  and  water  temperatures.  This  may  be 
explained  by  the  fact  that  increasing  the  air  and  water  tempera¬ 
tures  at  the  entrance  of  the  evaporation  tower  and  the  conden¬ 
sation  one  produces  more  thermal  energy  inside  the  solar  desa¬ 
lination  unit  which  in  its  part  increases  the  quantity  of  evaporated 
water.  This  leads  to  increasing  the  fresh  water  flow  rate. 

For  a  given  water  temperature  at  the  entrance  of  the  evapora¬ 
tion  tower  (Tc),  and  knowing  the  experimental  values  of  mc,  WC1 
and  Wc2,  the  humid  air  mass  velocity  is  [9]: 

(25) 

Using  this  equation,  the  air  circulates  in  the  unit  with  a  speed 
lower  than  the  theoretical  one.  Tests  show  that  GeXpS0.60  x  Gth- 
The  decrease  in  the  air  flow  rate  is  compensated  by  an  increase  in 
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Fig.  11.  Impact  of  solar  radiation  on  the  solar  collector  outlet  temperature. 


Fig.  13.  Cooling  water  temperature  variation  at  the  condensation  tower  exit  as  a 
function  of  water  temperature  at  the  evaporation  tower  inlet. 


Fig.  12.  Air  humidity  variation  at  the  top  of  the  evaporation  tower  as  a  function  of 
water  temperature  at  the  evaporation  tower  inlet. 


Fig.  14.  Variation  of  the  condensation  flow  rate  as  a  function  of  water  temperature 
at  the  evaporation  tower  inlet. 
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Fig.  15.  Variation  of  the  condensation  flow  rate  as  a  function  of  air  temperature  at 
the  condensation  tower  inlet. 


Table  3 

The  accuracy  of  the  experimental  results. 


Parameter 

Average 
relative  error 

ear  (%) 

Maximum 
absolute  error 
emax  (°C) 

Minimum 
absolute  error 
emin  (°C) 

Two 

2.468 

1.99 

3.5  x  i(r2 

Wg 2 

13.304 

16.36 

9.1  x  io-’ 

Te  2 

7.683 

4.28 

9.5  x  io-’ 

mc  (as  a  function  of  7i2) 

7.786 

5.31 

6.3  x  io-’ 

mc  (as  a  function  of  Tc 2 

i  4.347 

1.5 

4.7  x  itr2 

With:fmax  =  max  | Texp(i)~ 

^sim(0l>  emin  =  Him  |T 

exp(0— ^sim(0l>  I== 

:1,  2 

the  difference  between  the  air  temperatures  at  the  top  and  that  at 
the  bottom  of  the  condensation  tower.  This  is  why  the  actual  water 
distillation  rate  is  close  to  or  sometimes  higher  than  the  theore¬ 
tical  rate,  as  seen  in  Figs.  14  and  15. 

6.3.  Experimental  error 

Table  3  presents  a  simple  statistical  analysis  of  the  relative 
error  between  simulation  and  experimental  measurements.  The 
accuracy  of  the  experimental  results  may  be  calculated  using  the 
following  definition  [26]: 

100  *,  |Texp(i)-Tsim(i)| 

£ar(%)=  T  Texp(i)  (26) 

where  TeXp  denotes  the  experimental  data  of  T,  while,  Tsim  and  k 
are  the  theoretical  prediction  of  T  and  the  number  of  experimental 
measurements,  respectively.  The  error  analysis  of  the  different 
validated  parameters  shown  in  Figs.  11-15  is  presented  in  Table  3. 
It  is  clear  from  this  table  that  the  agreement  between  experiment 
data  and  the  results  calculated  from  the  theoretical  prediction  of 
the  present  device  is  fairly  good.  The  maximum  average  relative 
error  was  found  to  be  about  13%  for  the  absolute  air  humidty.  This 
value  was  expected  due  to  uncertainties  in  the  measured  quan¬ 
tities  of  the  experiments  and  the  simplifying  assumptions  in  the 
mathematical  models  in  addition  to  the  computational  errors  and 
the  accuracy  of  calculations. 

7.  Conclusion 

Brackish  water  or  seawater  desalination  can  be  an  interesting 
alternative  to  guarantee  a  secure  source  of  water.  Solar  energy  is 


considered  as  one  of  the  most  sustainable  way  to  supply  the 
energy  needs  for  desalination,  since  it  can  be  offered  close  to  the 
desalination  plants  where  connection  to  the  public  electric  grid  is 
either  not  cost  effective  or  not  feasible  and  avoid  environmental 
problems  associated  with  fossil  fuels.  Therefore,  the  improvement 
of  small  desalination  units  presents  an  attractive  solution  to 
decentralised  water  demand,  which  is  the  case  of  remote  and 
rural  area  in  developing  countries.  The  presented  unit  can  be 
considered  as  one  of  the  best  choice  for  production  of  fresh  water 
due  to  its  flexibility,  simplified  design,  low  maintenance,  extended 
life  time  for  over  twenty  years,  quasi  null  energy  consumption, 
construction  and  adaptation  for  use  in  rural  areas. 

In  this  paper,  a  steady  state  mathematical  model  expressing  the 
heat  and  mass  transfers  in  each  component  of  the  unit  is 
formulated.  The  developed  model  is  simulated  to  study  the 
influence  of  operating  conditions  on  the  thermal  performance  of 
the  unit.  It  was  found  that: 

•  The  GOR  increases  with  increasing  the  water  mass  flow  rate  at 
the  inlet  of  the  evaporation  tower.  There  is  an  optimum  value  of 
the  water  mass  flow  rate  where  the  GOR  is  maximum.  The 
optimal  water  mass  flow  rate  is  about  0.4  kg/s. 

•  The  GOR  decreases  with  increasing  water  mass  flow  rate  inside 
the  condensation  tower  and  then  stabilizes. 

•  The  thermal  efficiency  of  the  water  solar  collector  increases 
with  the  increase  of  the  water  mass  flow  rate. 

•  The  thermal  efficiency  of  the  water  solar  collector  decreases 
with  an  increase  of  the  inlet  water  temperature. 

•  The  thermal  efficiency  decreases  with  a  decrease  of  solar 
radiation  intensity. 

•  For  an  inlet  water  temperature  equal  to  the  ambient  tempera¬ 
ture,  the  solar  radiation  intensity  has  no  effect  on  the  thermal 
efficiency  of  the  water  solar  collector. 

Experimental  results  were  compared  with  the  simulation 
results.  It  was  shown  that  the  developed  models  are  able  to 
predict  accurately  the  trends  of  the  heat  and  mass  characteristics 
of  the  evaporation  and  condensation  chambers  and  the  solar 
collectors.  As  a  result,  the  proposed  models  can  be  used  for  sizing 
and  testing  the  behaviour  of  such  a  type  of  desalination  unit. 
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